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Abstract 

A new  so-called  truncation  error  reduction  method  (TERM)  is  developed  in 
this  work.  This  is  an  iterative  process  which  uses  a coarse  grid  (2h)  to  estimate 
the  truncation  error  and  then  reduces  the  error  on  the  original  grid  (h).  The 
purpose  is  to  use  multigrid  and  simple  stencils  for  high-order  accuracy. 

1.  Introduction 

The  driven  force  for  this  work  is  to  develop  a new  method  which  can  use 
simple  stencils  with  multigrid  method  for  high  order  accuracy.  The  multigrid 
(Brandt,  1984)  was  originally  used  to  accelerate  the  convergence  for  elliptic 
systems.  The  work  here  is  to  use  multigrid  for  high-order  accuracy  with  a simple 
stencil  which  allows  to  achieve  high-order  accuracy  with  much  fewer  points  than 
the  traditional  finite  schemes.  The  high-order  scheme  is  particularly  important 
for  direct  numerical  simulation  (Orszag,  S.  A.  & Patterson,  G.  S.,  1972;  Moine 
and  Mahesh,  1998)  and  the  large  eddy  simulation  (LES)  (Lilly,  1966;  Leseiur  & 
Metais,  1996). 

Let  us  take  a look  at  the  problem  with  coarse  grid  DNS  from  the 
mathematical  point  of  view.  Here,  the  coarse  grid  DNS  means  to  use  a grid  which 
is  acceptable  by  currently  available  computers.  The  problem  with  the  coarse  grid 
DNS  is  really  caused  by  the  truncation  error  while  the  mesh  size  ‘h’  and  time 
step  ‘k’  are  not  small  enough.  There  are  resolution  problems  with  coarse  grid 
DNS  as  well,  of  course.  Before  we  can  reduce  the  truncation  error,  we  need  to 
give  an  estimation  which  is  given  by  an  iterative  process  including  a coarse  grid 
discretization  and  a coarse-to-fine  grid  interpolation. 

We  use  the  Poisson  equation  as  a test  case  to  check  the  TERM  method  and 
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will  use  this  method  for  CFD  late. 

2.  Truncation  errors  reduction  methods  for  Poisson  equation 

Let  us  take  a look  at  the  following  Poisson  equation: 

+ = —k27t2(sin(k7tx)  + cos(lc7ty))  (1) 

dx  dy 

u(0,  y)  = cos(kJCy)',  m(1,  y)  = cos(k7ty) ; (2) 

u(x,0 ) = 1 + sin(&;cc);  m(jc,1)  = (-1)*  + sin  (km) 

The  exact  solution  is  sin(km)  + cos(&7Ty) . We  adopt  a standard  second 


order  central  difference  scheme  with  uniform  grids,  Ax  = Ay  = h,  and  obtain  a 
finite  difference  scheme: 


u — 2m  4-  u 1 u — 2m  1 4*  u 

Ui-Uj  zui,j+uM,j  , ui,H  ^ui.j+uU+ 1 , „h 

.2  .2  Si,j 


h2 


h2 


where  u is  the  exact  solution  and  s(  . = —k 2tz 2 sin(A'7tri  )~k27t 2 cos {k7Cxi ) . 

This  finite  difference  scheme  will  give  an  exact  solution  if  we  can  find  the 

exact  truncation  error,  T*  • , which  has  a second  order,  0(h2) . We  all  know  it  is 

difficult  to  find  the  exact  truncation  error,  but  we  can  try  to  give  a more  accurate 
estimate  by  an  iterative  process  which  can  be  described  as  follows. 

1 .  Assume  T*’°  = 0 , do  one  or  two  times  Gauss-Seidel  iterations,  multigrid 


(Brandt,  1984)  may  be  much  faster,  to  get  an  approximation u-  'j . 

2.  Find  the  approximation  on  the  coarse  grid  uj 2hj2  = /^w^, . In  our  case,  a 

simple  injection  is  used  for  if/'  (or  uf£j2  = unji  )• 

3.  Estimate  the  truncation  error  on  the  coarse  grid: 
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.2/1,1 


= s 


t2h 


-2  uv]+u2h 


>.j 


M ,/ 


.2/1 


(2/I)2 


, - 2u  +uin 

i,j~i  ni,j  i,j 

(2hf 


(4) 


4.  Estimate  the  truncation  error  on  the  fine  grid: 


,/i,i 


1 


T — _Ly  t 
S'l.yl  “ .*l2hLi2,j2 


(5) 


where  /*,,  should  have  same  or  higher  order  than  the  finite  difference 


scheme  and  we  should  use  1/16  instead  of  1/4  in  the  formula  if  we  use  a 
fourth-order  scheme. 

5.  Get  a revised  finite  difference  scheme  on  the  fine  grid: 


4-,-2^+C.  . 

p + p + <« 

6.  Go  to  step  1 and  start  a new  loop  until  II  xh '"+1  -xh'n  II  < tolerance 


This  method  can  be  extended  with  some  changes  to  other  flow  problems  which 
are  usually  non-linear  and  time-dependent. 


3.  Numerical  test  for  Poisson  equation  with  different  grids 

A number  of  numerical  tests  by  using  this  truncation  error  reduction  method 
(TERM)  for  different  grids  (8x8,  16x16,  32x32,  64x64,  128x128)  and  different  k 
( 1 , 2,  3,  4,  8,  1 6)  have  been  conducted. 

Table  1 compares  the  numerical  results  with  the  exact  solution.  The 
comparison  shows  the  error  ratio  of  coarse  grid  2h  over  fine  grid  h for  standard 


second  order  central  difference  scheme. 


2h 

/ 

„ h 

u—u 

/ 

u — u 

, is  near  4,  which 


shows  the  numerical  scheme  has  a second  order  accuracy  (Burden  & Faires, 
1996).  By  using  the  TERM  method  , the  accuracy  has  been  significant  increased. 
First  we  find  the  error  ratio  is  nearly  16,  which  shows  we  got  a fourth  order 
accuracy  with  second  order  stencils.  From  Table  1 , we  can  find  a solution  with  a 
grid  of  16x16  by  the  new  method  can  get  same  results  as  one  obtained  with  a grid 
of  64x64  by  the  standard  central  difference  scheme.  It  shows  we  can  save  the 
grids  or  the  computer  memory  by  16  times  for  this  2-d  Poisson  equation.  The 
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computation  time  can  be  saved  almost  100  times.  Similar  conclusion  can  be 
found  with  other  grids  and  different  wave  number  k.  More  encouraging,  the 
method  shows  a same  achievement  when  the  wave  number  becomes  higher.  It 
still  can  improve  the  results  when  we  have  only  three  points  for  one  wave  ( see 
the  cases  of  8x8  for  k=4  and  16X16  for  k=8).  It  provides  a high  potential  which 
could  be  used  for  coarse  grid  DNS  for  a more  accurate  numerical  simulation  for 
transitional  turbulent  flow. 

Figure  1-4  show  the  L2  norm  of  errors  between  the  numerical  solution  and 

exact  solution  against  the  standard  central  difference  method  and  the  TERM 
method  for  different  grids.  It  clearly  shows  that  the  new  method  increases  the 
accuracy  significantly  and  has  potential  for  more  accurate  coarse  grid  simulation. 
For  2-D,  the  new  method  saves  at  least  16  times  in  grid  point  numbers  and  much 
more  in  CPU  time.  We  can  anticipate  it  will  save  at  least  64  times  in  the  number 
of  grid  points  for  3-D  problems. 

4.  Conclusion 

The  new  so-called  truncation  error  reduction  method  (TERM)  can 
significant  increase  the  accuracy  of  numerical  solution  for  Poisson  equation  for 
different  wave  numbers.  It  shows  TERM  can  use  coarse  grids  to  achieve  much 
more  accurate  numerical  solution  than  the  standard  finite  difference  scheme  for 
Poisson  equations.  It  can  save  hundreds  times  in  memory  and  computational  time. 
Potentially,  the  method  may  be  used  for  fluid  dynamics  and  for  coarse  grid  DNS. 
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Analytic  Solution 
u{x,  y)  = sink7Dc  + cos  kjcy 

IMAX  = 8 
JMAX  = 8 

IMAX  = 16 
JMAX  = 16 

IMAX=32 

JMAX=32 

IMAX=64 

JMAX=64 

IMAX=I28 

JMAX=128 

3.8297E-003 

2.6273E-004 

1 .6809E-005 

7.8555E-004 

6.1889E-005 

4.3033E-006 

2.8284E-007 

1.81 16E-008 

aMaaasMWi 

3.578  IE-002 

9.3823E-003 

2.4108E-003 

6.1159E-004 

1.5406E-004 

9.9472E-003 

8.4777E-004 

6.0961  E-005 

4.0573E-006 

2.61 10E-007 

9.3301  E-002 

2.3893E-002 

6.1033E-003 

1.546  IE-003 

3.8931  E-004 

ID 

4.2942E-002 

3.8647E-003 

2.8834E-004 

1.944  IE-005 

1.2565E-006 

1.9922E-00I 

4.8970E-002 

1.2387E-002 

3.1303E-003 

7.8776E-004 

BiMM 

1.1337E-002 

8.7461  E-004 

5.9764E-005 

3.8789E-006 

2.34 19E  001 

5.5260E-002 

I.3733E-002 

3.4418E-03 

ID 

6.7739 

1.6488E-001 

1.2195E-002 

8.9128E-004 

5.9201  E-05 

■■■jaiwav 

wSjSSgrnswwiW 

91.4400 

1.5034 

2.5141E-001 

5.8164E-002 

1.4329E-02 

ERROR 

(TERM) 

92.3567 

31.4476 

I.8054E-00I 

1.268  IE-002 

9.0120F.-04 

homogeneous  isotropic  turbulence,  Phys.  Rev.  Lett.  Vol  2,  pp76-79,  1972. 


Table  1 . Error  comparison  between  central  difference  scheme  and  TERM 
with  different  grids  and  wave  numbers 


JMAX 

Figure  I : Comparison  of  errors  between 
standard  Central  difference  scheme  and 
TERM  scheme,  k=l . 


■MAX 

Figure  2:  Comparison  of  errors  between 
standard  Central  difference  scheme  and 
TERM  scheme,  k=2. 
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Figure  3:  Comparison  of  errors  between 
standard  Central  difference  scheme  and 
TERM  scheme,  k=3. 


Figure  4:  Comparison  of  errors  between 
standard  Central  difference  scheme  and 
TERM  scheme,  k=4. 


Appendix 

We  can  prove  the  above  TERM  method  for  Poisson  equation  has  fourth 
order  for  interior  points.  For  simplicity,  we  consider  the  one  dimensional  case: 


a 2U  = 

dx2  J 


(7) 


Similarly  to  (3),  the  standard  second  order  central  difference  scheme  with 
uniform  grids,  Ax  = h , can  be  written  as 


M,_l  —2m,.  +um  _ 

,.2  "*’ti  Ji 


(8) 


♦ injecting  points 

fine  grid:  • ♦ • ♦ • ♦ • 

• interpolating  points 


coarse  grid:  • 


Figure  5 Injecting  points  and  interpolating  points 


1.  Injecting  points  where  the  fine  grid  point  coincides  with  the  coarse  grid  point 

(u2h 

On  coarse  grid 
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Ui-2  - 2lii  + ui+2  , _ _ , 

(2ft)2  Wi 


(9) 


Here,  T(-  is  the  truncation  error  on  coarse  grid.  We  have  used  the  following 


relation  between T,  andT( 


ifr  Ui-2~2Ui  +«/+2  V » 


T,  = — T-  + f.  = 

' 4 ' '4 


4/r 


+ T,.  (10) 


7 


Here,  f(.  is  the  truncation  error  for  the  new  scheme.  Substitute  (10)  into  (8),  we 
have 


— w(._  2 +16m,_,  — 30w;  +16miX|  — m, 


*+l  , S' 


12/t2 


+ r,=f,  0D 


Using  Taylor  series,  we  can  find  that  T,  ~ 0(h  ‘ ) and  (11)  is  a fourth  order 
scheme. 

2.  Interpolating  points  where  the  fine  grid  point  is  located  between  two  coarse 
grid  points. 

On  coarse  grid 


(2/t)2 


(12) 


2l/[  + | +M,  + 3 _ 

(2/t)2  i+‘  Jm 


(13) 


Here,  T.  , and  Ti+1  are  the  truncation  errors  on  coarse  grid.  We  have  used  the 


following  relation  forT,_,  ,T(+I  and  T. 


1 1 


T.  = - X — X (T  + T . ) + f . = - 

* 4 2 8 


/,_i  + /,+!  - 


^1-3  Ui- 1 *0+1  + *0  + 3 

Ah2 


+ T, 
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(14) 

Here,  T.  is  the  truncation  error  for  the  new  scheme.  Substitute  (16)  into  (8),  we 
have 

— m,_3  + 33mm  —64 Uj  + 33m, +l  — uj+3  ^ (/,■_ i ~2/.  + f i+l)  | ~ _ y 
24/t2  6 


We  have  following  relations  on  the  fine  grid: 


h 2 

m,_ | - 2m,.  + Mi+1  h 


-* (4,-,  )=/,-, 

(16) 

II 

Mjp: 

i 

_i 

(17) 

(18) 

, and  x,.  < £,+1  < x/+2 . 

Substitute  (16), 

h 2 


ui  - 2m, +l  + m;^, 

/t2 


(17),  (18)  into  (15),  we  have 

— m,_3  +4m,_2  +17mm  -40m,-  +17m,+1  + 4m,+2  -m,.+3  | ~ 

24h2  ' 

1^-  («<4>  ) - 2«  <->  («, ) + «<4>  («,  .,»  = /, 


(19) 


Because  m<4) (£_,  ) - 2m(4) (£ ) + m(4)  (&+i  ) « m(6)  (r?,  )/t2 


(20) 


here  x,._2  < t],  < x,+2 . Therefore 

^ (m<4> (£M ) - 2m<4> (!,.)  + » t4) (£,., ))  - 0(/!4)  (21) 

Using  Taylor  series,  we  can  find  that  T,  ~ 0(h4)  and  (15)  is  a fourth  order 
scheme. 


